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(3) A method and apparatus for calibrating magnetic and geometrical distortions in an imaging system. 

(57) A method and apparatus to calibrate and 
correct magnetic and geometrical distortions in 
an imaging system, specifically a computer to- 
mography (CT) system, includes a calibration 
object which is placed on the surface of an 
image intensifies An image frame of the calib- 
ration object is generated, and a vertical correc- 
tion table is generated corresponding to the 
amount of image distortion caused to the image 
frame of the calibration object The image frame 
is corrected using the vertical correction table, 
and the corrected image frame data is summed 
or averaged and either stored for further trans- 
formation into a slice image or used to generate 
a horizontal correction table to correct distor- 
tions in the horizontal direction. Once the verti- 
cal and horizontal correction tables have been 
generated using the calibration object, subse- 
quent image frames are corrected using the 
correction values in the correction tables. 
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The present invention relates to a method and apparatus for calibrating and correcting magnetic and 
geometrical distortions in an imaging system, and, more specifically, to calibrating and correcting magnetic and 
geometrical distortions in an image intensrfier-based computed tomography (CT) system. 

5 Description of the Related Art 

An image intensifier-based CT system utilizes an image intensifier to detect x-rays emitted by an x-ray 
source after they are passed through an object under test. The image intensifier converts x-ray photons to opti- 
cal photons to electrons, which are accelerated and focused down the length of the image intensifier onto a 
w small phosphor screen. The phosphor screen generates an image which is recorded by an optical sensor which 
outputs a video signal typically formatted into 525 lines of video information. This video signal is then digitized 
for further processing. 

A two-dimensional slice image of an object under test is reconstructed from a plurality of the image frames 
corresponding to a plurality of views of the object under test taken as either the object under test is rotated 
15 360° on a turntable (turntable rotation) or as the x-ray source and image intensifier are rotated 360° around 
the stationary object under test (gantry rotation). For example, if the system produces 30 image frames per 
second and the turntable or gantry completes one rotation every 60 seconds, the slice image will be recon- 
structed from 1800 separate image frames corresponding to the 1800 distinct views of the object under test. 
Image intensifiers impart two distortions to the positions of detected x-ray photons. The first distortion is 
20 due to the construction geometry of the image intensifier and is referred to as "geometrical" distortion. This dis- 
tortion causes a pagoda-like distortion of the detected image frame. That is, if a set of nested squares are formed 
by wires and placed on the entrance surface of the image intensifier, the output image appears as shown in 
Fig. 1 . This is also referred to as "pin-cushion" distortion. 

The second distortion is caused by the magnetic force of the earth, and will be found in any imaging system 
25 using an x-ray detector having an electron tube. Because the x-ray detector (image intensifier) accelerates elec- 
trons down its length, these electrons are subject to magnetic distorting forces of the form: 



30 



F = q v x B f 



(1) 



where F is the distorting force vector; 

q is the elementary charge of the electron; 

— • 

v is the velocity vector of the electron; and 

B is the magnetic field vector of the earth. This distorting force further twists the image of Fig. 1. 

If the imaging system is a turntable rotation system, these two distortions are constant for each image frame 
generated by the image intensifier. However, if the imaging system is a gantry rotation system, the magnetic 
distortion differs for each image frame generated by the image intensifier. The magnetic distorting force distorts 
each image frame differently due to the different orientation of the image intensifier as the gantry rotates. 

Additional magnetic distortions result in any imaging system using a tube TV camera as an optical sensor 
as a result of the moving electrons in the TV imaging tube. 

SUMMARY OF THE INVENTION 

Accordingly, it is an object of the present invention to provide a method and apparatus for calibrating and 
correcting magnetic and geometrical distortions in an imaging system. 

It is another object of the present invention to provide a method and apparatus to calibrate and correct ver- 
tical magnetic and geometrical distortions in an imaging system. 

It is another object of the present invention to provide a method and apparatus to calibrate and correct hori- 
zontal magnetic and geometrical distortions in an imaging system. 

It is another object of the present invention to provide a method and apparatus to calibrate and correct ver- 
tical and horizontal magnetic and geometrical distortions in an imaging system. 

It is another object of the present invention to provide a calibration object which can be used to calibrate 
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imaging systems to correct for magnetic and geometrical distortions. 

Additional objects and advantages of the invention will be set forth in the description which follows, and in 
part will be apparent from the description, or may be learned by practice of the invention. The objects and advan- 
tages of the invention may be realized and obtained by means of the instrumentalities and combinations par- 
5 ticularly pointed out in the appended claims. 

To achieve the foregoing objects, and in accordance with the purposes of the invention as embodied and 
broadly described herein, an apparatus is provided for correcting distortions in an imaging system, comprising: 
a calibration object having at least one reference line; a device for generating an image frame of the calibration 
object; a device for generating a correction table including a plurality of correction values corresponding to 
w amounts of distortion of the reference line at a plurality of locations in the image frame; and a device for cor- 
recting distortions of an object under test using the plurality of correction values in the correction table. 

To further achieve the foregoing objects, and in accordance with the purposes of the invention as embodied 
and broadly described herein, a method is provided for correcting distortions in an imaging system using a cali- 
bration object having at least one reference line, comprising the steps of: generating an image frame of the 
15 calibration object; generating a correction table including a plurality of correction values corresponding to 
amounts of distortion of the reference line at a plurality of locations in the image frame; and correcting distortions 
of an object under test using the plurality of correction values in the correction table. 

BRIEF DESCRIPTION OF THE DRAWINGS 

20 

The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate 
a presently preferred embodiment and method of the invention and, together with the general description given 
above and the detailed description of the preferred embodiment given below, serve to explain the principles of 
the invention. Of the drawings: 
25 Fig. 1 is a drawing illustrating "pin-cushion" distortion in an image-intensifier based CT system; 

Fig. 2 is a block diagram of a distortion correcting imaging system in accordance with the present invention; 

Fig. 3 is a diagram of a calibration object in accordance with the present invention; 

Fig. 4 illustrates a displayed image frame of the calibration object of Fig. 3, which has not been corrected 

for distortions; 

30 Fig. 5 is an expanded view of 32lines (a partial frame) of the imageframe in Fig. 4, with only one horizontal 

line being illustrated; 

Fig. 6 is a graph showing a parabolic interpolation to detect a minimum data value; 
Fig. 7A is a vertical correction table before being normalized in accordance with the present invention; 
Fig. 7B is a vertical correction table after being normalized in accordance with the present invention; 
35 Fig. 8A illustrates which data values are corrected in accordance with the vertical distortion correction tech- 

nique of the present invention; 

Fig. 8B illustrates a displayed image of the vertically corrected partial image frame of Fig. 5; 
Fig. 9 is a preliminary table utilized for the horizontal correction technique of the present invention; 
Fig. 1 0A is a horizontal correction table before completion in accordance with the present invention; and 
40 Fig. 10B is a horizontal correction table after completion in accordance with the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

Reference will now be made in detail to a presently preferred embodiment of the invention as illustrated in 
45 the accompanying drawings. 

Fig. 2 illustrates a CT imaging system in accordance with the present invention. The CT imaging system 
includes x-ray source 10, x-ray detector 12, optical sensor 14, video selector 16, A-to-D converter 18, buffer 
20, CPU 22, and display 24. Optionally, the system also includes collimator 26. Although buffer 20 is shown 
as a separate component, more typically buffer 20 constitutes part of the CPU 22. 
so An object under test 28 is disposed on a turntable (not shown) for rotation between the x-ray source 10 

and the x-ray detector 12 in accordance with conventional techniques. Alternatively, the x-ray source 10 and 
the x-ray detector 12 are disposed on a rotating gantry (not shown) which is rotated around the stationary object 
under test 28 in accordance with conventional techniques. 

X-ray source 10 emits a cone of x-rays toward the object under test 28. The x-rays may be collimated by 
55 collimator 26 to select the slice width of the object under test 28, or the entire cone of x-rays may be permitted 
to pass through the object under test 28 to the x-ray detector 12. 

As embodied herein, x-ray detector 12 comprises an image in tensifier which, as discussed, converts x-ray 
photons to optical photons to electrons, which are focused and accelerated toward a phosphor screen 30 con- 
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tained in the image intensifier 12. The phosphor screen 30 reconverts the electrons back to optical photons to 
form an image frame of the object under test 28. Preferably this image frame is 525 lines of video information. 

The image frame on the phosphor screen 26 is recorded by optical sensor 14. As embodied herein, optical 
sensor 14 may comprise a camera, a charge-coupled diode converter or some other solid state device suitable 
5 for recording the image frame from the phosphor screen 30 of the image intensifier 12. In a preferred embodi- 
ment, the optical sensor 14 comprises an MTV-32D saticon T.V. camera, manufactured by Toshiba Corp. 

Video selector 16 is provided to select pixels, lines and frames of video information output by the optical 
sensor 14 for digitization in A-to-D converter 18. Typically, an A-to-D converter digitizes each of the 525 video 
lines into 512 pixels of video information. Copending U.S. Patent Application No. 593,261 filed on even date 
10 herewith and assigned to the same assignee as this application (the entire contents of which are hereby incor- 
porated by reference), discloses a preferred embodiment of video selector 16. As described therein, video 
selector 16 can be programmed through CPU 22 to select the total number of frames of video information to 
be digitized, the specific frames out of the total number of frames generated during a scanning period to be 
converted, the specific lines in each frame to be converted, and the number of pixels in each line to be con- 
15 verted. This selected video information is converted to digital information by the digital converter 18. 

As is also described in the above mentioned copending patent application, the video selector 16 also pro- 
vides a control signal to the x-ray source 10 to provide dose control overthe amount of x-rays to which an object 
under test is exposed. This is particularly important if the object under test is a human being or an animal. 

The selected digitized data is transferred to the buffer memory 20 for temporary storage before being for- 
20 warded to the CPU system memory (not shown) in the CPU 22. CPU 22 processes the digitized data to correct 
magnetic and geometrical distortions and to provide CT slice images which are displayed on display 24. 

CPU 22 may comprise any analog or digital computer or computational device or devices having sufficient 
memory allocation and operating speed to carry out the distortion correction techniques of the present invention. 
In a preferred embodiment CPU 22 comprises a personal computer, such as the AST-386 PC, which may be 
25 modified to provide additional processing power as necessary. 

As embodied herein, display 24 is provided with inputs from both the optical sensor 14 and the CPU 22. 
In this implementation, a switching device (not shown) is provided to determine which source of video input 
will be displayed. Alternatively, separate displays could be provided for separately displaying the output video 
information of the optical sensor 14 and the CPU 22. 
30 The correction techniques of the present invention provide for distortion corrections in both the vertical and 

horizontal directions. 

In accordance with the calibration and correction techniques of the present invention, before generating 
image frames of the object under test, image frames of a calibration object are generated to "measure" both 
the vertical and horizontal distortions. 

35 As shown in Fig. 3, a calibration object 32 comprises a grid of opaque objects 34. As embodied herein, 

the opaque objects comprise copper wires 0.5 mm in diameter. A first set of 27 parallel wires, 0.5 inches apart, 
are provided in one direction, and a second set of 27 parallel wires, 0.5 inches apart are provided perpendicular 
to the first set of 27 parallel wires. The wires are contained within a block of acrylic plastic 36, 14" on each side 
and 11/16" thick. The dimensions of the wire grid, wires and plastic block should be chosen to meet and accom- 

40 modate the dimensions of the image intensifier 12. 

Both vertical and horizontal corrections are performed in CT imaging systems using either a rotating table 
configuration or a rotating gantry configuration. Corrections for systems using gantry rotation are more complex 
than for systems using rotating turntables because, as discussed, each image frame taken as the gantry rotates 
around the object under test exhibits different distortions due to the magnetic forces operating on the system, 

45 which change with the gantry rotation. Vertical and horizontal correction techniques of the present invention 
will first be described for a CT system using a rotating table to rotate the object under test with respect to a 
stationary x-ray source and x-ray detector. These correction techniques are henceforth referred to as "turntable 
corrections." 

so 1. Turntable corrections 

Vertical distortion correction is provided by placing the calibration object 32 on the face of image intensifier 
12 and generating an image frame of the calibration object 32. Alternatively, the calibration object can be dis- 
posed on a turntable (not shown). However, the translation technique utilized to generate the slice image would 
55 have to be accommodated for the spacing between the turntable and the image intensifier 12, in accordance 
with conventional techniques. 

A display of the image frame of the calibration object might look as shown, for example, in Fig. 4. Due to 
geometrical and magnetic distortions, the horizontal grid wires 40 of the calibration object are not perfectly hori- 
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zontal. Further, the vertical wires 42 in the center of the image frame appear closer together than those at the 
edges of the image frame. The vertical correction technique corrects the bending of the horizontal grid wires 
and the horizontal correction technique corrects the spacing between the vertical grid wires. 

Through experimentation, it has been Found that if 32 lines of the 525 lines comprising the image frame 

5 of the calibration object are selected and displayed, the complete distortion of at least one horizontal grid wire 
of the calibration object 32 can be observed. If these 32 lines are selected from the central portion of the image 
frame, it can be assumed for purposes of this invention, that the horizontal grid wires are distorted equally by 
the distorting forces and appear parallel (although bent) in this central portion of the image frame. Thus, col- 
limator 26 and video selector 16 are adjusted to select a predetermined set of at least 32 lines from the central 

10 portion of the image frame to implement the vertical correction technique. 

To more clearly illustrate the vertical distortion correction, Fig. 5 shows an expanded displayed image of 
32 lines selected from the central portion of the image frame of the calibration object 32. For convenience, these 
32 lines have been renumbered 0-31. Henceforth, these 32 lines will be referred to as a "partial image frame." 
Also for convenience, only one horizontal grid wire 40 is illustrated in Fig. 5. 

15 If this image of the calibration object 32 were not distorted, the horizontal grid wire 40 would appear as a 

straight line across, for example, line 12 of the displayed partial image frame. However, due to distortions, as 
illustrated, the horizontal grid wire 40 is distorted by eight lines at the left edge of the partial image frame (column 
0). Further, for each column of the partial image frame, there exists an offset value corresponding to the number 
of lines by which the horizontal grid wire 40 is displaced from line 12. For example, at columns 0-3 the offset 

20 from line 12 is eight lines (20-12), whereas at column 256, the offset from line 12 is four lines (16-12). This 
offset is calculated for each column to generate a vertical correction table, and the data values of each column 
of subsequent image frames or partial image frames are offset by a correspond ing amount to correct for vertical 
distortions in the image frame or partial image frame. 

To determine the offset for each column, the location of the distorted horizontal wire 40 (in terms of image 

25 frame lines) is identified for each column. As embodied herein, this is accomplished by determining the average 
data value of all 32 pixels in each column. Data values corresponding to the grid wires will typically be signifi- 
cantly lower than the average data value because the grid wires block x-rays from being detected at that loca- 
tion. To prevent identifying erroneous data values as corresponding to the horizontal grid wire 40, as embodied 
herein, a predetermined threshold is selected (for example one-half the dynamic range of the data values in 

30 the column, below the column average) below which data values are known to correspond to the grid wires. 
Unless the data value is below this threshold, the system will not consider it to correspond to horizontal wire 
40. Thus, the CPU begins at column 0, line 0 to determine if the data value for that pixel is a predetermined 
threshold level below the average data value for column 0. If not, the CPU reads the data value at line 1 of 
column 0 to make the same determination, and so on until the horizontal grid wire 40 is identified in column 0. 

35 As embodied in Fig. 5, when the CPU reads the data value at line 20 of column 0, it identifies this as being 

the horizontal grid wire 40. The CPU then compares the data value of line 20 to that of line 21 to determine 
which data value is lower so as to determine the center of the horizontal grid wire 40. For example, although 
in Fig. 5, the horizontal grid wire 40 is shown as being one-dimensional, in reality the horizontal grid wire 40 
has some width and may comprise data values of varying intensities from, for example, lines 19-21. The CPU 

40 continues to compare subsequent data values to locate more precisely the line corresponding to the center of 
the horizontal grid wire 40. 

In a preferred embodiment, a parabolic interpolation is further employed to identify the theoretical center 
of the horizontal grid wire 40 which may occur between image frame lines. For example, as shown in Fig. 6, a 
graph of the discrete data values for lines 0-31 of column 0 corresponds to the solid curve 60. However, as 

45 shown by the dotted curve 62, a parabolic interpolation of these data values may identify the theoretical mini- 
mum of the parabola (hence, the center of the grid wire) at line 20.3. This decimal precision can be used to 
provide a weighting function to grade subsequent corrections, as will be described below. Any conventional 
method or technique for providing this interpolation may be employed in accordance with the invention. 

Thus, as embodied herein, the CPU identifies the horizontal grid wire 40 as corresponding to line 20.3 in 

so column 0. This value is stored in a memory location of a vertical correction table corresponding to column 0. 

This same procedure is repeated 51 1 times, once for each additional column of data values in the partial 
image frame, until a vertical correction table of 512 values has been generated. This table is shown, for example, 
in Fig. 7A, which has not been normalized. 

The table of values of Fig. 7A is then normalized to zero as shown in Fig. 7B, to generate a vertical correction 

55 table having offset values for each column of the partial image frame. This correction table can be used to cor- 
rect vertical distortions for any subsequent image frame or partial image frame generated by the CT image sys- 
tem. For example, having generated the vertical correction table shown in Fig. 7B, an operator may desire to 
provide a slice image of an object under test. Assuming that the system generates 30 image frames each sec- 
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ond and that the rotating table supporting the object under test makes one complete rotation in 60 seconds, 
1800 image frames will comprise the CT slice image. The vertical correction table of Fig. 7B is utilized to correct 
the data values for each of these 1800 image frames or partial image frames. 

Specifically, the image frame lines are corrected using the offset values of the vertical correction table of 
5 Fig. 7B in accordance with the following relationship: 

DoutM] = (1-F).D in [c,f + "I + F D in M + 1+1] (2) 

where D 0Ut Ic,f] = the corrected data value for the pixel corresponding to column (c) and line (?); 

where D ln [c,£] = the uncorrected data value for the pixel corresponding to column (c) and line (f); 

10 where I = the integer portion of the correction value for column (c) taken from the vertical cor- 

rection table; and 

where F = the fractional portion of the correction value for column (c) taken from the vertical cor- 

rection table. 

As shown in Fig. 8A, correctionof the partial image frame of Fig. 5, in accordance with relationship (2) 

is above, has the effect of selecting only data values from the cross-hatched portion of the uncorrected partial 
image frame for the vertically corrected data set. 

-The vertical correction technique of the invention reduces the number of data values in the vertically cor- 
rected data set. This is illustrated, for example in Fig. 8B, which shows a display of the corrected partial image 
frame of Fig. 8A. As shown, only lines 0-23 (the cross-hatched portion) remain in the data set corresponding 

20 to the corrected partial image frame. As illustrated, the bend in the horizontal grid wire 40 is straightened out 
in accordance with the correction values of the vertical correction table. 

As discussed above, the fractional portion F of the offset va lue s in the ve rtical correction table are utilized 
in the correction technique of the present invention to provide a weighting function between data values of adja- 
cent lines of an image frame or partial image frame. For example, in accordance with relationship (2) above, 

25 the corrected data value corresponding to column 0, line 0 of the corrected partial image frame shown in Fig. 
8B, is .7 times the data value from column 0, line 8 of the distorted partial image frame shown in Fig. 8A. plus .3 
times the data value for column 0, line 9 of the partial image frame shown in Fig. 8A. In this manner, the fractional 
portions of the correction values in the vertical correction table are utilized to grade the intensity of the corrected 
data values based on adjacent line data values. 

30 Alternatively, a faster computational implementation is possible with a nearest neighbor correction tech- 

nique rather than the linear interpolation technique provided in relationship (2) above. In the nearest neighbor 
correction technique, the table correction values are rounded to the nearest whole number and the vertical cor- 
rection is provided in accordance with the relationship: 

DooJc,*] = D ln [c,£ + Q, (3) 

35 

where D^Icf ] = the corrected data value for the pixel corresponding to column (c) and line (I); 

where D| n [c,f ] = the uncorrected data value for the pixel corresponding to column (c) and line (2 ); and 

where I = the nearest whole number to the correction value for column (c) taken from the vertical 

correction table. 

40 The 24 corrected partial image frame lines are then summed or averaged to provide a single line of 512 

data values. 

This process is repeated 1800 times, once for each of the 1800 image frames which comprise the slice 
image, to generate 1800 single lines of 512 data values. Preferably, this process is performed in real time to 
reduce the overall data load. 

45 These data values may be displayed as what is commonly known as a sinogram. However, more typically, 

these data values are further processed and manipulated in accordance with conventional CT techniques to 
generate the slice image of the object under test. As embodied herein, the inverse Radon transform implemen- 
ted using the convolution backprojection algorithm is utilized to transform the data values to slice image form. 
However, any conventional transformation techniques may be utilized to transform the data to slice image form 

so for display. 

Although the slice image may be comprised of only vertically corrected image frames or partial image 
frames, more typically, the vertically corrected image frames or partial image frames are further corrected for 
horizontal distortions before the slice image is displayed. Horizontal correction of the distorted image frames 
or partial image frames will now be discussed with respect to the turntable corrections. 
55 As embodied herein, the horizontal correction technique of the present invention corrects the perceived 

irregular spacing between the vertical grid wires 42, which are distorted from their regular interval in the calib- 
ration object 32. The horizontal correction is performed on the sinogram data values after the image frames or 
partial image frames have been vertically corrected and summed or averaged into single lines of data values 
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for each image frame or partial image frame of the slice image. However, a horizontal correction table must 
first be generated. 

For horizontal correction, data values corresponding to the distorted vertical grid wires 42 of the calibration 
object 32 must be identified. As discussed, these data values are lower than the average data value of a pre- 
5 determined image frame line of data values, and a predetermined threshold is selected to dearly identify the 
vertical grid wires and prevent identifying spurious noise changes as incorrect vertical grid wire locations. 
The system scans a line of the vertically corrected and summed or averaged data values to detect the col- 
__ umn locations corresponding to vertical grid wires 42 of the calibration object 32. This is done in a manner similar 
to the detection of horizontal grid wires 40 in each column for the vertical correction technique. However, as 
10 embodied herein, the scan is made along only one line of the sinogram data values for the horizontal correction 
technique. All 512 columns are scanned to identify the locations of each vertical grid wire 42. Also in a manner 
similar to that performed for the vertical correction technique, a parabolic or other type of interpolation technique 
my be employed to identity the theoretical centers of the vertical grid wires 42. 

Since the calibration object 32 embodied herein has 27 vertical grid wires, the system will identify 27 col- 
ts umns corresponding to the vertical grid wire centers. These 27 columns are stored in a preliminary table. This 
is illustrated in Fig. 9, where, for example, the data value for the first vertical grid wire in the calibration object 
has been found to occur at column 4.2, the data value for the second vertical grid wire has been found to occur 
at column 25.7, up through the data value for the 27th vertical grid wire which has been found to occur at column 
510.1. 

20 Using the 27 column location values from the preliminary table, a horizontal correction table is generated. 

There are 26 intervals between the first and 27th vertical grid wires. Dividing the 512 pixel locations of a line 
of the sinogram by 26 intervals yields a value of 1 9.69 as the number of pixel locations or columns per interval. 
As embodied herein, for computational simplicity and efficiency, only the integer portion of this value (19) is 
used. Multiplying this value (1 9) times the number of intervals (26) yields 494 as the number of correction values 

25 in the horizontal correction table embodied herein. Thus, one value from the preliminary table of 27 values is 
located every 19 locations in the horizontal correction table, as shown in Fig. 10A. Thus, in the preferred embo- 
diment, the horizontal correction table has a lower number of values (494) than the number of columns (512). 

Linear interpolation is then used to fill in the intermediate values of the horizontal correction table according 
to the relationship: 

30 i 

F(19n + i) = • (fn + 1 " ^n) + fn (4) 



where F (19n+i) = correction table values between the given 27 values from the preliminary table; 
where i = the location in the horizontal correction table between two given values from the prelimi- 

35 nary table 0^i^1 8; and 

where f n = . given 27 values in the preliminary table 0^n^25. 

For example, the interpolated value for column two in the horizontal correction table is calculated as: 



40 



F2 « • (25.7 - 4.2) + 4.2 

19 



= 6.5 

45 

This process is repeated until the horizontal correction table is completed, as shown for example in Fig. 
10B. 

Each line of the vertically corrected and summed or averaged data values is then horizontally corrected 
using the correction values from the horizontal correction table in accordance with the following relationship: 
50 D out [c] = (1-F) • D jn [l] + F - D in [l + 1] (5) 

where D out [c] = the corrected data value corresponding to column (c); 
where D in [l] = the uncorrected data value corresponding to column (I); 

where I = the integer portion of the correction value for column (c); and 

55 where F = the fractional portion of the correction value for column (c). 

Only data values occurring between the columns corresponding to vertical grid wires 1 and 27 are considered 
for the horizontal correction. 

Similar to the vertical correction, the fractional portions of the correction values in the horizontal correction 
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table are utilized to grade the intensity of the corrected data values based on adjacent column data values. 

Alternatively, the nearest neighbor correction technique can be used to provide the horizontal correction. 
In this embodiment, the table values are rounded to the nearest whole number, such that relationship (5) is 
reduced to: 

5 D^c] = D ln fl] (6) 

where D out [c] = the corrected data value corresponding to column (c); 
where D ln [l] = the uncorrected data value corresponding to column (I); and 

where I = the nearest whole number to the correction value for column (c). 

10 As each line of data values corresponding to the 1 800 image frames (the vertically corrected sinogram infor- 

mat ion) are horizontally co rr ected, the y are storedjo generate a vertically and horizontally corrected set of data 

values corresponding to sinogram information. These data values are then further processed and transformed 
using the inverse Radon transform implemented by the convolution backprojection algorithm into a slice image 
as described above for the vertical correction technique of the present invention. The slice image is then dis- 

15 played on display 24. 

2. Gantry Corrections 

For imaging systems utilizing a gantry to rotate the x-ray source 10 and the x-ray detector 14 with respect 
20 to a stationary object under test 28, the vertical and horizontal corrections are similar to the turntable correc- 
tions. However, separate vertical and horizontal correction tables must be generated for each of the 1 800 image 
frames which comprise the slice image. This is because, as the gantry rotates through 360° around the object 
under test, at each gantry position corresponding to each image frame, the magnetic distorting forces acting 
on the system differ. 

25 For example, for the preferred embodiment in which the gantry makes a complete rotation around the object 

under test every 60 seconds and takes 30 image frames per second, totaling 1 800 image frames per scan, the 
vertical correction table for gantry correction requires 1800 separate lines of correction values, one for each 
image frame. Thus, the vertical correction table becomes two-dimensional, having 1800 lines and 512 correc- 
tion values per line. Each line of the two-dimensional vertical correction table is generated in a manner similar 

30 to that described above with respect to the turntable correction. The process is simply repeated 1 800 times, 
once for each image frame. 

In an alternative embodiment the vertical correction table may comprise only a fractional portion of the 
number of frames collected during a scan, such as 225, 450 or 900 separate lines of correction values. Each 
line is used to correct a corresponding number of image frames, such as 8, 4 or 2, respectively. 

35 After the vertical correction table has been generated, subsequent distorted image frames are corrected 

in a manner similar to the vertical turntable correction described above. However, each image frame is corrected 
by the corresponding line of the two-dimensional vertical correction table. 

In the alternative embodiment in which a lesser number of lines are used in the vertical correction table, 
more than one image frame or partial image frame is corrected using the same line of the vertical correction 

40 table. For example, rf the vertical correction table is 225 lines, the first 8 image frames of the gantry rotation 
are corrected using line 1 of the vertical correction table, the second 8 image frames are corrected using line 
2 of the vertical correction table, and so on. 

As with the vertical turntable correction technique, the corrected image frame lines for each image frame 
are summed or averaged into a single line which comprises the sinogram data values to be transformed into 

45 the slice image. These vertically corrected and summed or averaged data values are used to perform the hori- 
zontal gantry correction. 

To perform the horizontal gantry correction, a horizontal correction table is generated in a manner similar 
to that described with respect to the horizontal turntable correction technique. However, because of the gantry 
rotation, the horizontal correction table is again two-dimensional and preferably has one line of correction values 
so for each image frame, i.e. 1800 lines of 494 correction values. Each line of the two-dimensional horizontal cor- 
rection table is generated in a manner similar to that described above with respect to generating the horizontal 
correction table for horizontal turntable corrections. The process is simply repeated 1 800 times, once for each 
image frame. 

As with the vertical correction table, the number of correction lines in the horizontal correction table can 
55 be varied such that one line is used to correct more than one image frame. 

After the horizontal correction table has been generated, subsequent distorted image frames are corrected 
in a manner similar to the horizontal turntable correction technique. However, for horizontal gantry corrections, 
each image frame is corrected by the corresponding line of the horizontal correction table. If the horizontal cor- 
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rection table includes less than one line for each image frame, more than one image frame is corrected using 
the same line of the horizontal correction table. 

The data values, vertically and horizontally corrected, comprise sinogram data values which are transfor- 
med into the slice image in accordance with conventional techniques, as described above with respect to turnt- 
5 able corrections. 

In the foregoing embodiments corresponding to both turntable and gantry corrections, vertical and hori- 
zontal corrections are performed. However, in an alternative embodiment, it is not necessary to perform the 
.vertical correction, and only horizontal correction is performed. In particular, collimator 26 can be used to select 
the width of the x-ray fan beam to be narrow enough such that detected x-rays corresponding to the distorted 

10 x-ray fan beam lie within the selected number of frame lines selected by the video selector. In this embodiment, 
if the data values for each column of the resulting image frame or partial image frame are summed or averaged 
to a single line of video data values, the image frame need not be vertically corrected. However, a horizontal 
distortion correction in accordance with the present invention may still be performed. 

While there has been illustrated and described what are at present considered to be preferred embodiments 

15 and methods of the present invention, it will be understood by those skilled in the art that various changes and 
modifications may be made, and equivalents may be substituted for elements thereof without departing from 
the true scope of the invention. In addition, many modifications may be made to adapt a particular element, 
technique or implementation to the teachings of the present invention without departing from the central scope 
of the invention. Therefore, it is intended that this invention not be limited to the particular embodiments and 

20 methods disclosed herein, but that the invention include ail embodiments and methods falling within the scope 
of the appended claims. 



Claims 

A method of correcting distortions in an imaging system using a calibration object having at least one refer- 
ence line, comprising the steps of: 

generating an image frame of said calibration object; 

generating a correction table including a plurality of correction values corresponding to amounts of 
distortion of the reference line at a plurality of locations in said image frame; and 

correcting distortions of an object under test using said correction values in said correction table. 

A method according to claim 1, wherein said step of generating an image frame of said calibration object 
includes the steps of: 

irradiating said calibration object with x-rays; 
detecting x-rays which pass through said calibration object; 
generating a set of data values corresponding to said detected x-rays; and 
selecting a predetermined set of said data values corresponding to a predetermined slice width. 

A method of correcting distortions according to claim 1, or claim 2 wherein said image frame of said cali- 
bration object includes a plurality of columns and lines of data values, and the step of generating a cor- 
rection table includes the steps of: 

identifying locations of the reference line in each column of said image frame; 
storing a plurality of values corresponding to the identified locations of the reference line for each 
column of said image frame; and 

normalizing to zero the stored plurality of values corresponding to the identified locations of the 
reference line. 

A method according to claim 3, wherein the step of identifying includes the step of detecting a data value 
below a predetermined minimum. 

A method according to claim 4, wherein the step of identifying further includes the step of interpolating 
between data values. 

A method according to claim 1, wherein the step of correcting the distortion of an object under test com- 
prises the steps of: 

generating an image frame of said object under test, said image frame including a plurality of col- 
umns and lines of data values; and 
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offsetting data values of each column of said image frame of the object under test by an amount 
equal to a corresponding correction value in the correction table. 

7. A method according to claim 6, wherein the step of offsetting is performed for each column of said image 
5 frame of the object under test in accordance with the relationship: 

DautM] - (l-FJ-DmM + I] + F D ln [c,f + 1+1] 

where DouJcC] = the corrected data value corresponding to column (c) and line (2 ); 

where Djn[c t £] = the uncorrected data value corresponding to column (c) and line (2 ); 

10 where I = the integer portion of the correction value for column (c) taken from the vertical 

correction table; and 

where F = the fractional portion of the correction value for column (c) taken from the ver- 

tical correction table. 

15 8. A method according to claim 6, wherein the step of offsetting is performed for each column of said image 
frame of the object under test in accordance with the relationship: 
DoJcf] = D ln [c,f + I], 

where D out [c,f ] = the corrected data value for the pixel corresponding to column (c) and line (2 ); 

20 where D m [c,2] = the uncorrected data value for the pixel corresponding to column (c) and line 

(0; 

where I = the nearest whole number to the correction value for column (c). 

9. An apparatus for correcting distortions in an imaging system, comprising: 
25 a calibration object having at least one reference line; 

means for generating an image frame of said calibration object; 

means for generating a correction table including a plurality of correction values corresponding to 
amounts of distortion of the reference line at a plurality of locations in said image frame; and 

means for correcting distortions of an object under test using said correction values in said correc- 
30 tion table. 

10. An apparatus according to claim 9, wherein said calibration object comprises a grid of opaque lines having 
a first set of parallel lines in a first direction and a second set of parallel lines perpendicular to said first 
set of parallel lines. 



35 



11. An apparatus according to claim 10, wherein said opaque lines of said calibration object comprise wires. 



12. An apparatus according to claim 9, wherein said means for generating an image frame of said calibration 
object includes: 

40 means for irradiating said calibration object with x-rays; 

means for detecting x-rays which pass through said calibration object; 

means for generating a set of data values corresponding to said detected x-rays; and 

means for selecting a predetermined set of said data values corresponding to a predetermined slice 

width; 

45 

13. An apparatus according to claim 9, wherein said image frame of said calibration object includes a plurality 
of columns and lines of data values, and said means for generating a correction table includes: 

means for identifying locations of the reference line in each column of said image frame; 
means for storing a plurality of values corresponding to the identified locations of the reference line 
so for each column of said image frame; and 

means for normalizing to zero the stored plurality of values corresponding to the identified locations 
of the reference line. 

14. An apparatus according to claim 13, wherein the identifying means includes means for detecting data 
55 values below a predetermined minimum. 

1 5. An apparatus according to claim 1 4, wherein the identifying means further includes means for interpolating 
between data values. 

10 
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16. An apparatus according to claim 9, wherein said means for correcting the distortion of an object under 
test includes: 

means for generating an image frame of said object under test, said image frame including a 
plurality of columns and lines of data values; and 
5 means for offsetting data values of each column of said image frame of the object under test by an 

amount equal to a corresponding correction value in the correction table. 

17. An apparatus according to claim 16, wherein the offsetting means includes means for offsetting the data 
values for each column of said image frame of the object under test in accordance with the relationship: 

w D out [c,f] = (1-F)D ln [c,f + I] + F D in M + 1+1] 



where D olrt [c,f ] = the corrected data value corresponding to column (c) and line (I); 

where Dm[c,P] = the uncorrected data value corresponding to column (c) and line (?); 

where I = the integer portion of the correction value for column (c) taken from the vertical 

is correction table; and 

where F = the fractional portion of the correction value for column (c) taken from the ver- 

tical correction table. 

18. An apparatus according to claim 16, wherein the offsetting means includes means for offsetting the data 
20 values for each column of said image frame of the object under test in accordance with the relationship: 

D ou t[c,£] = D in [c.£ + I], 



where DouJc.f ] = the corrected data value for the pixel corresponding to column (c) and line (2 ); 

where D^c,?] = the uncorrected data value for the pixel corresponding to column (c) and line 

25 (f): 

where I = the nearest whole number to the correction value for column (c). 

19. A method of correcting distortions in an imaging system using a calibration object having a plurality of sub- 
stantially parallel reference lines, comprising the steps of: 
30 generating an image frame of said calibration object; 

generating a correction table including a plurality of correction values corresponding to amounts of 
distortion between adjacent reference lines of the calibration object; and 

correcting distortions of an object under test using said correction values in said correction table. 

35 20. A method according to claim 1 9, wherein said step of generating an image frame of said calibration object 
includes the steps of: 

irradiating said calibration object with x-rays; 

detecting x-rays which pass through said calibration object; 

generating a set of data values corresponding to the detected x-rays; and 
40 selecting a predetermined set of said data values corresponding to a predetermined slice width. 

21. A method according to claim 19, wherein said image frame of said calibration object includes a plurality 
of columns and lines of data values, and said step of generating a correction table includes the steps of: 

identifying the columns of the data values corresponding to the reference lines; 
45 storing a column number corresponding to each identified reference line at equally spaced locations 

of a correction table; and 

filling-in the locations of the correction table between the column numbers. 

22. A method according to claim 21 , wherein the step of filling-in includes the step of linearly interpolating be- 
50 tween the stored column numbers: 

23. A method according to claim 19, wherein the step of correcting the distortion of an object under test 
includes the steps of: 

generating an image frame of said object under test, said image frame including a plurality of col- 
55 umns and lines of data values; 

summing or averaging the data values of each column of said image frame of said object under 
test into a single data value for each column so as to generate a single line of summed or averaged data 
values; and 
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offsetting the data values of selected columns of the image frame of the object under test by an 
amount equal to a corresponding correction value in the correction table. 

24. A method according to claim 23, wherein the step of offsetting is performed in accordance with the rela- 
5 tionship: 

OoAc] = <1-F)Din[l] + F.D ln [l + 1] 

where D^c] = the corrected data value corresponding to column (c). 

where D^l] = the uncorrected data value corresponding to column (I) 

10 where I = the integer portion of the correction value from the horizontal correction table cor- 

responding to column (c); and 

where F = the fractional portion of the correction value from the horizontal correction table cor- 

responding column (c). 

15 25. A method according to claim 23, wherein the step of offsetting is performed in accordance with the rela- 
tionship: 

Doutlc] = DM 

where D^c] = the corrected data value corresponding to column (c); 

20 where D^fl] = the uncorrected data value corresponding to column (I); and 

where I = the nearest whole number to the correction value for column (c). 



26. An apparatus for correcting distortions in an imaging system, comprising: 

a calibration object having a plurality of substantially parallel reference lines; 
25 means for generating an image frame said calibration object; 

means for generating a correction table including a plurality of correction values corresponding to 
amounts of distortion between adjacent reference lines of the calibration object; and 

means for correcting distortions of an object under test using said correction values and said cor- 
rection table. 

30 

27. An apparatus according to claim 26, wherein said means for generating an image frame of said calibration 
object includes: 

means for irradiating said calibration object with x-rays; 
means for detecting x-rays which pass through said calibration object; 
35 means for generating a set of data values corresponding to the detected x-rays; and 

means for selecting a predetermined set of said data values corresponding to a predetermined slice 

width. 

28. An apparatus according to claim 26, wherein said image frame of said calibration object includes a plurality 
40 of columns and lines of data values, and further comprising: 

means for summing or averaging the data values of each column of said image frame of said cali- 
bration object into a single data value so as to generate a single line of summed or averaged data values; 
and 

wherein said means for generating a correction table includes: 
45 means for identifying the columns of the data values corresponding to the reference lines; 

means for storing a column number corresponding to each identified reference line at equally 
spaced locations of a correction table; and 

means for filling-in the locations of the correction table between the column numbers. 

50 29. An apparatus according to claim 28, wherein said means for filling-in includes means for linearly interpolat- 
ing between the stored column numbers: 

30. An apparatus according to claim 26, wherein said means for correcting the distortion of an object under 
test includes: 

55 means for generating an image frame of said object under test, said image frame including a 

plurality of columns and lines of data values; 

means for summing or averaging the data values of each column of said image frame into a single 
data value for each column so as to generate a single line of summed or averaged data values; 

12 
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means for offsetting the data values of selected columns of the image frame of the object under 
test by an amount defined by a corresponding correction value in the correction table. 

31. An apparatus according to claim 30, wherein the means for offsetting the data values includes means for 
5 offsetting the data values in accordance with the relationship: 

DoJc] = (1-F) . DM + F • D ln p + 1] 

where D ou t[c] = the corrected data value corresponding to column (c). 

where D|n[l] = the uncorrected data value corresponding to column (I) 

10 where I = the integer portion of the correction value from the horizontal correction table cor- 

responding to column (c); and 
"where F = thefractional portion of the correction value from the horizontal correction table cor- 

responding to column (c). 

15 32. An apparatus according to claim 30, wherein the means for offsetting the data values includes means for 
offsetting the data values in accordance with the relationship: 

DoJc] = DJI] 

where D out [c] = the corrected data value corresponding to column (c); 

20 where D h [l] = the uncorrected data value corresponding to column (I); and 

where I = the nearest whole number to the correction value for column (c). 



33. A method of correcting distortions in an imaging system using a calibration object having at least one refer- 
ence line in a first direction and a plurality of substantially parallel reference lines in a second direction, 

25 comprising the steps of: 

generating an image frame of said calibration object; 

generating a first correction table including a plurality of correction values corresponding to amounts 
of distortion of the reference line at a plurality of locations in said image frame; 

correcting distortions of an object under test using said correction values in said first correction ta- 

30 ble; 

generating a second correction table including correction values corresponding to amounts of dis- 
tortion between the plurality of substantially parallel reference lines; and 

correcting distortions of the object under test using said correction values in said second correction 

table. 

35 

34. An apparatus for correcting distortions in an imaging system, comprising: 

a calibration object having at least one reference line in a first direction and a plurality of substan- 
tially parallel reference lines in a second direction; 

means for generating an image frame of said calibration object; 
40 means for generating a first correction table including a plurality of correction values corresponding 

to amounts of distortion of the reference line at a plurality of locations in said image frame; 

means for correcting distortions of an object under test using said correction values in said first cor- 
rection table; 

means for generating a second correction table including correction values corresponding to 
45 amounts of distortion between the plurality of substantially parallel reference lines; and 

means for correcting distortions of the object under test using said correction values in said second 
correction table. 

35. A method of correcting distortions in an imaging system using a calibration object having at least one refer- 
so ence line, said imaging system including a rotating gantry, said method comprising the steps of: 

rotating a gantry around said calibration object; 

generating a plurality of image frames of said calibration object as said gantry rotates around said 
calibration object; 

generating a correction table having a plurality of lines of correction values, each line of correction 
55 values corresponding to ampounts of distortion of the reference line at a plurality of locations in one of 

said plurality of image frames; and 

correcting the distortions of a plurality of image frames of an object under test using said correction 
values from said correction table, wherein one of said plurality of lines from the correction table is used 

13 
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to correct at least one of said plurality of image frames of said object under test. 

36. An apparatus for correction distortions in an imaging system including a rotating gantry, said apparatus 
comprising: 

5 a calibration object having at least one reference line; 

means for rotating said gantry about said calibration object; 

means for generating a plurality of image frames of said calibration object as said gantry rotates 
about said calibration object; 

means for generating a correction table having a plurality of lines of correction values, each line of 
10 correction values corresponding to amounts of distortion of the reference line at a plurality of locations in 

one of said plurality of image frames; and 

means for correcting the distortions of a plurality of image frames of an object under test using said 
correction values from said correction table, wherein one of said plurality of lines from the correction table 
is used to correct at least one of said plurality of image frames of said object under test. 

15 

37. A method of correcting distortions in an imaging system using a calibration object having a plurality of sub- 
stantially parallel reference lines, said imaging system including a rotating gantry, said method comprising 
the steps of: 

rotating a gantry around said calibration object; 
20 generating a plurality of image frames of said calibration object as said gantry rotates around said 

calibration object; 

generating a correction table having a plurality of lines of correction values, each line of correction 
values corresponding to amounts of distortion between adjacent reference lines of the calibration object 
for one of said plurality of image frames; and 
25 correcting distortions of a plurality of image frames of an object under test using said correction 

values from said , correction table, wherein one of said plurality of lines from said correction table is used 
to correct at least one of said plurality of image frames of said object under test. 

38. An apparatus for correcting distortions in an imaging system including a rotating gantry, said apparatus 
30 comprising: 

a calibration object having a plurality of substantially parallel reference lines; 
means for rotating said gantry around said calibration object; 

means for generating a plurality of image frames of said calibration object as said gantry rotates 
around said calibration object; 
35 means for generating a correction table having a plurality of lines of correction values, each line of 

correction values corresponding to amounts of distortion between adjacent reference lines of the calib- 
ration object for one of said plurality of image frames; and 

means of correcting distortions of a plurality of image frames of an object under test using said cor- 
rection values from said correction table, wherein one of said plurality of lines from said correction table 
40 is used to correct at least one of said plurality of image frames of said object under test. 
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